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(<~^ . We consider a model where, for a finite disintegrating system, yields of 
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composites can be calculated to arbitrary accuracy. An analytic answer for 
yields is also known in the thermodynamic limit. In the range of tempera- 



> 

. ture and density considered in this work, the model has a phase transition. 



This phase transition is first order. The analytic expression for yields of 
composites, in the thermodynamic limit, does not conform to the expression 
< ria >= a^'^ f [a" {T — Tc)) where, the usual identification would be that Tc is 

b . 

^ . the critical temperature and r, a are critical exponents. Nonetheless, for finite 

systems, we try to fit the yields with the above expression. A minimisation 
procedure is adopted to get the parameters Tc, r and a. While deviations 
from the formula are not negligible, one might believe that the deviations are 
consistent with the corrections attributable to finite particle number effects 
and might then conclude that one has deduced at least approximately the 
values of critical parameters. This exercise thus points to difficulties of trying 
to extract critical parameters from data on nuclear disintegration. An inter- 
esting result is that the value of Tc deduced from the "best" fit is very close to 
the temperature at which the first order phase transition occurs in the model. 
The yields calculated in this model can also be fitted quite well by a 
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parametrisation derived from a droplet model. 
25.70.Pq, 24.10.Pa, 64.60.My 
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I. INTRODUCTION 



The following parametrisation, often used to fit nuclear multifragmentation data, gives 
an elegant expression for yields of composites: 

<na>=a--f{a'^{T-T,)) (1.1) 

Here a is the mass number of the composite, Tc is the critical temperature, r is the Fisher 
exponent and the expansion is valid in the neighbourhood of Tc and for "large" a. Variants 
of the equation are also used. For example, in the percolation model one replaces T — by 
p — Pc where pc is the value of p at which an infinite cluster first appears . The values of 
T and a for three dimensional percolation model are 2.18 and 0.45 respectively. In the case 
of the lattice gas model, also used in a great deal to study nuclear multifragmentation 0, 
Eq (1.1) is valid not only in the neighbourhood of Tc but along the entire Kertesz line and 
extends also to lower than critical density 

There is also another class of models which have been used to fit multifragmentation 
data. Typical of these is the Copenhagen statistical multifragmentation model 0. Details 
vary between different versions but the principal assumption is that, at some larger than 
normal volume, hot nuclear matter breaks up into different clusters. Nuclear interactions 
between different clusters are deemed insignificant. The coulomb interaction between differ- 
ent composites can be taken into account at least approximately. Numerous applications of 
such models have been made with impressive success. The advantage of such models is that 
quantum effects such as shell effects can be included. 

How well does Eq (1.1) work for such models? This is a relevant question in view of the 
fact that there already exist many applications of such models to fit actual data f^^- We 
take a simplified version of such models so that exact answers can be obtained (as opposed 
to Monte-Carlo sampling which make calculations very long). The story that unfolds is quite 
interesting and is the subject of this paper. For simplicity, we first consider a model of one 
kind of particle in the next section. This is then anlysed for parametrisation. The method 
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of analysis is explained in section 3. Results are presented in section 4. After looking at fits 
to models with isotopic spins with and without the Coulomb force, we return to the model 
of one kind of particles again, this time, trying a fit with a droplet model. A comparison 
between different models is presented in section 7. Our conclusions are summarised in section 
8. 



II. MODEL OF ONE KIND OF PARTICLE 



Details of the model can be found in [T^JTll • They are summarised here for completeness. 
We first consider just one kind of particles as the thermodynamic properties for this model 



were studied in detail ||10| , p!2| -|l^ . We can have monomers or composites of k nucleons. The 
composites have ground state energy —Wk + a(T)k'^^^. The first term is the volume energy 
with 14^=16 MeV. The second term is the surface tension term. The surface tension term is 
taken to be temperature dependent as in 0: a{T) = ao[{T^ -T^)/{T^ + T^)]^/^ with ao =18 
MeV and = 18 MeV. The internal partition function of a composite of A; > 1 nucleons is: 

Zk = exp[(Wk - a{T)e/^ + T^k/eo)/T] (2.1) 

where we have used the standard Fermi-gas model for excited states(the third term on the 
right hand side of Eq.(2.1)). For A; = 1 we take z\=\. 

We want to construct the canonical partition function of the system which has total A 
nucleons in a volume V which is much larger than the volume of a nucleus of A nucleons in 
the ground state:\^ > Vq. We assume when clusters are formed in the volume V ^ they do 
not overlap with each other. Thus the available volume free to the clusters is less than V . 
It is given by V/ = l^ — Vexc- We take Vexc to be Vq = A x 6 where h = l/.16/m~'^. In reality 
Vexc should also depend on multiplicity but this complication is ignored here. The canonical 
partition function Qa_{T) of A nucleons is then given by 

Qa = T.U^ (2-2) 

fc>i 

where LOk is the partition function of one composite of k nucleons: 



ujk = ^{2TxmTf/^k'/^ X Zk (2.3) 



and the sum rule must be obeyed:^ = Y^kuk- As noted before ||l0[, the partition function 
Qa for A nucleons can be easily generated in the computer by utilising a recursion relation. 
Starting with Qq = 1 one can build all higher ones using 

1 ^ 

Qp = - V kukQp-k (2.4) 
The expression for the yield of composites is, of course, of primary interest. This is given by 

< Uk >= u;fe— — (2.5) 
Several things are known for this model. The critical temperature for the model is the 



temperature at which the surface tension vanishes [0 and hence it is at T = Tc = 18MeV. 



The critical volume is at = Vq. At temperatures below 18 Mev there is first order 



phase transition [TyjT^. The temperature of phase transition depends upon the density. 
This temperature was called boiling temperature in [|10|. In the temperature range we are 
concerned with in this paper, there is only first order phase transition. The phase transition 
temperature is characterised by a sudden jump in specific heat. In finite systems we will 
take boiling temperature to be the temperature at which the specific heat maximises. 

We will try to fit the yields of Eq.(2.5) by the generic formula (Eq.(l.l)). The exact 
expressions (Eqs.(2.2) to (2.5)) give no clue of a simple parametrisation. Provided V is 
reasonably bigger than Vq (see however |T^) and A is also large we can use the grand 
canonical results to obtain some insight on simple parametrisation. This answer is well- 
known: 

< nj, >= exp[(/i(T) + W + T^le^)k/T - a{T)k^'^ /T] (2.6) 

There is no exact correspondence between Eq.(2.6) and Eq.(l.l). Thus we may at best hope 
to get an approximate fit. How we do it is detailed in the next section. 
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III. THE FITTING PROCEDURE 



Here we follow very elegant methods given in 0. For later use in the text, we will give 
adequate details. 

The quality of fit is given by the smallness of x^. If the calculated quantity F is a function 
of two parameters Y = Y{ai,bi) and we are trying to fit it to a function g{ai,bi,a, (3,'j...) 
then = J2iO^{0'i, h) — gicii, h, a, P, 7----))^/ '^i^j^ ^i)^- Variations of this criterion are 
also possible. Eq. (1.1) requires us to find "best" values for r, cr and Tc. This is done in 
several steps. 

(1) At each temperature we find the "best" r associated with an attempted fit < Ua >= 
a~'^C where C is a constant. This follows from Eq.(l.l) but only at T = T^, hence one can 
argue that if Eq.(l.l) were exact for yields calculated by Eq.(2.5) we would get null and 
the correct r at T = T^. This would determine both and r. Of course, null is not 
found since "exact" fit is not given by Eq.(l.l) However we can draw a best "r" vs. T curve 
for a pure power law. This r as a function of T will have a minimum which we label Tmin- 
Since at T^, the fit is strictly a power law, one can accept that temperature as where the 

of the fit is minimum. However, we will determine T^. using the method described in step 

3. 

(2) Let z = a'^{T — T^); f{z) has a maximum at some value of z = z: fmax = f{z)- For 
each mass number a the yield < na > as a function of temperature has a maximum at some 
value of temperature, Tmax{,o). At this temperature < na> {max) = a~'^ fmax where fmax is 
a constant independent of a. This relationship allows us to choose "best" values for r and 

fmax • 

(3) The value of r found from step (2) is higher than Tmin found at step (1). This means 
that if we look for T appropriate for r, two values of T are available from the graph at step 
(1). The lower value is to be chosen as the value of Tc. The scaling property (see steps (4) 
and (5)) is badly violated for the other choice. 

(4) Now that we know and Tmaxi^), the temperaure at which the yield of composite 
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a is maximised, we find, by least squares fit the "best" value of a from the condition 
0-'^ {Tmax{0') — Tc) = const. for all a. 

(5) The scaling law can now be tested by plotting < Ua > vs. a'^{T — Tc). Plots for 
all a's should fall on the same graph. 

In following steps (1) to (5) the range of a is to be chosen judiciously. It can not be very 
small (since Eq.(l.l) applies to "large" a's). But a should also be truncated on the high 
side (significantly smaller than the size of the dissociating system). 

IV. RESULTS FOR ONE KIND OF PARTICLES 

We present our results in Figs.l to 3. The sizes of the systems are taken to be A = 174 
and A = 240. The upper pannels of the figures use the freeze-out volume V = SVq and the 
lower pannels use V = AVq. Both are shown here for completeness. It will suffice here to 
discuss only the cases with V = SVq. 

Fig. 1 shows r vs. T drawn according to step (1) of section 3. The dotted line is the 
value of T deduced from step (2). This cuts the curve of step (1) at two temperatures (step 
(3)). The lower value of the temperature is taken as T^. In the same Fig. we also plot C^/A 
as a function of T. The peak of this curve corresponds to the first order phase transition, 
called "boiling" temperature in [0. It is remarkable that the of Eq.(l.l) is very close 
to "boiling temperature" in this particular example. In this respect, the model is similar to 
the behaviour of the lattice gas model where parametrisation Eq.(l.l) works on the lower 
side of the critical density also provided Tc is replaced by Tb, the temperature of first order 
phase transition ||^,|^ at the given density. In Fig.l we have also plotted the value of 
as a function T (step 1). In Fig.2 we plot In < Ua > vs. Ina. Two graphs are shown for 
each disintegrating systems. The graph with higher values of < > (shown as diamonds) 
follows from step (2) of section 3. These are the maximum values oi < Ha > for each a 
obtained at corressponding temperatures Tmax{o)- The lower values oi < Ha > (shown as 
stars) are all at the same temperature, namely at Tc which, for example, is 6.32 MeV for 
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A=17A and is 6.54 MeV for A=240. This is how r would be estimated from experiments 



T^. The crucial testing of the scaling law is shown in Fig.3. where we plot < Ha > 



vs. a"{T — Tc). For the range of a chosen (10 to 40) the results nearly fall on the same 
graph. Since one does not know a priori how much error is due to finite particle number of 
the disintegrating system, one might be tempted to to conclude that the fit to Eq. (1.1) is 
adequate. The parameters r, a from best fits are 2.72, 1.06 respectively for A=174 and 2.78, 
1.23 respectively for A=240. The deduced Tc are 6.32 MeV and 6.54 MeV which are very 
different from the critical temperature of 18 MeV for the model but compare remarkably 
well with the temperatures where the specific heats peak and which correspond to first order 
phase transition temperatures at the given densities. 

Here we want to comment that, as seen from Fig. 1, has a minimum at a temperature 
very close to Tc. So one may conclude that the methods of determining Tc using step 1 or 3 
yield almost the same result. 

These then are the two salient features: (1) Numerical fits of Eq.(l.l) are surprisingly 
close and (2) interpreting as the critical temperature is wrong although the deduced 
does correspond to a phase transition temperature. We briefly look at how these features 
hold for a more realistic model with isotopic spin. 



V. MODELS OF TWO KINDS OF PARTICLES 

We considered the following simplified version. Protons and neutrons are elementary 
particles. For deuteron, triton, ^He and "^He, experimental binding energies are used but 
no excited states are included. For masses greater than 4 a semi-empirical formula is used. 
Composites of n neutrons and z protons have binding energies W{n + z) — a{n + zY^^ — 
tts — o,z (n+ly/'-i • Here 1^=15.8 MeV, cr=18 MeV is the surface tension term (taken here 
as temperature independent), = 23.5MeV is the symmetry energy term and az=-72 MeV 
is the Coulomb energy term. For composites >4 we also include excited states in the Fermi- 
gas model (as in section II). We also incorporate Coulomb interaction between composites 
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in the Wigner-Seitz approximation 0. 

It is however also instructive to consider two kinds of particles with the surface tension 
and symmetry energy term but with the Coulomb energy switched off. Such matters, in 
the thermodynamic limit can be self-bound at zero temperature. With Coulomb included 
however matter would fall apart in the thermodynamic limit even at zero temperature. We 
show first the results for finite systems with surface and symmetry energy terms but no 
Coulomb. For brevity we will show results only for the choice V = SVq. For two kinds of 
particles, Eq.(l.l) can be used in more than one way. The subscript a in Eq.(l.l) can stand 
for the mass number a = n + z. Then Eq.(l.l) gives the distribution of particles of given 
mass number, irrespective of n or z. The lower panels of Fig. 4 to 6 use that. Or, as is more 
suitable for experiments, a can stand for z, the charge, irrespective of the neutron number n. 
The upper panels of Figs. 4 to 6 use that. We notice that the deduced value of continues 
to be close to where the specific heat maximises. 

For brevity, with Coulomb included, we show the results for V = 3Vo only. Distribution 
of particles with a given charge, irrespective of the neutron number n are shown in Figs. 7 
to 9. The most noticeable difference with the "no Coulomb" cases is that the deduced is 
now farther removed from the temperature at which the specific heat maximises. Coulomb 
energy introduces rather significant changes to the caloric curves. This was noticed before 
||TT| . For example, without the Coulomb, as the particle number increases, the peak of the 
specific heat becomes narrower and the temperature at which the maximum is obtained 
slightly shifts to higher value (see Fig. 4). The narrowing and the shifting of the peak seems 
to be almost exactly compensated by the long range repulsive Coulomb force as more and 



more charged particles are added (Fig. 4 here and also Fig. (6) in [|TT| where this is discussed 
in greater detail). Properties of nuclear matter but where protons carry charges, require 
further studies. 
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VI. FIT TO A DROPLET MODEL 



We will now try to fit the predicted yields given by the model of section II (one kind 
of particles) with a well-known droplet model An early application of the model, to 



heavy- ion collisions, can be found in [16|. The model has been revived recently ||17|| . 

In Fisher's model, the condensation of a real gas into large drops(clusters) of liquid is 
modeled. This shares various similar features with multifragmentation models. For instance, 
the potential energy of large clusters consists entirely of a bulk term and a term associated 
with the loss of binding energy at the surface. There is no Fermi energy term, as the 
molecules inside the cluster are assumed to be Boltzmann distributed. The entropy of large 
clusters is, however, more complicated. As clusters become large, the dominant effect may 
be ascribed, once again, to a bulk term, and the remainder to a surface term. It was 
pointed out that liquid clusters may not be restricted to spherical shapes, as is the case 
in most multifragmentation models. This prohibits the use of one form of surface area to 
parametrize the surface contribution. In reference it is argued that, at low temperatures, 
the most important configurations will be compact and globular. Their surface areas s are 
not much greater than the minimum possible, and are assumed to admit the asymptotic 
condition that, 

s{k,P)/k^O as k oo and 
s{k, (3)/ log k ^ oo as A; ^ oo. (6.1) 

Where, k is the number of molecules occupying the cluster. If, for finite clusters, one 
introduces this surface area to calculate various surface contributions, one must introduce 
a correction term which varies as rlogfc. The sign and magnitude of r is estimated from 
various other considerations involving other models. 

One may thus, very generally, obtain the mean number of clusters of size k as 

< Uk >= Ck-^ exp((/ig - fii)k/T + csfc'/Vr) (6.2) 
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Here both fig and fii are functions of T. At coexistence and also at the critical temperature, 
they become equal to one another. Also C2 is a function of temperature and at the 
coefficient C2 goes to zero. Since above Tc there is no distinction between the liquid and 
the gas phase, one can not speak of droplets. Thus the theory only applies to T < Tc. As 
such the formulation is more limited than the model of Eq. (1.1) which applies to both 
sides of Tc. The following fit was tried. We set r = 2. let a = (/ig — fJ'i)/T, 7 = C2/T. 
We fit the calculated < nfc > to Ck~'^ exp{ak + 7^^/^) at different temperatures where 
a, 7 values at each temperature are varied for best fit. The values of 0,7 as a function of 
temperature is shown in Fig. 10 where we also show rather remarkable fit with the values of 
< Uk > obtained from the model of section II. The values of a and 7 both go to zero near 
temperature T = 6.5 MeV suggesting that the critical temperature is 6.5 MeV. 



VII. RELATIONSHIP BETWEEN DIFFERENT MODELS 

In this section a grand canonical approximation to the model of section II is constructed. 
As in the previous section, this will be concentrated in the region below the boiling temper- 
ature. A different parametrization of the output from the exact canonical calculation will 
be obtained. Though this parametrization has not been as throughly investigated as that of 
the previous section, it will serve to provide a qualitative understanding of the behaviour of 
the yields below the boiling temperature. 

The yields(Fig. (|Tl|) , Fig.(0)) are obtained from the exact expressions Eq.( |2.5|) . For the 
present analysis they are adequately approximated by Eq.( p.6|) . The parametrization offered 
by Eq. (|2.6|) is of the form 

(rik) = Ck^'^ exp(aA; + 7^2/^) (7.1) 

The above is different from the parametrization of the droplet model where r = 2. Though 
the expressions look similar, the fit parameters C,a,7 will assume values different from 
those of C, a, 7. The interpretation of a is also quite precise in this approximation. 
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When using Eq. (|2.6|) to estimate the result of a canonical calculation, the free parameter 
is the chemical potential fi. It is usually determined by imposing that the model correctly 
reproduce the total number of particles composing the system(i.e., J2k{^k)k = A ). This 
is a complicated problem in general. A clue may be obtained by observing the behaviour 
of /i as obtained from the canonically calculated Helmholtz free energy F = —TlogQ^. A 
plot of /i obtained thus is plotted in Fig. (|l^) (data points). The behaviour of this /i may 
be estimated by the following simple argument. Far below the boiling temperature the 
system exists mostly as one large cluster and a few small ones(see Fig. (0)). The large 
cluster is considered as the liquid state. Far above the boiling temperature the system 
exists mostly as many small clusters: this is considered as the gas phase. The chemical 
potential of either system may be estimated by keeping the system in contact with a heat 
reservoir, adding one particle to the system, and noting the change in free energy, i.e.. 



/i 



AF 



= F{T,V,A + 1) — F{T,V,A). On entering the system, the new particle, 



may apriori attach itself to any of the existing clusters, or simply thermalize as a monomer 
in the system. It will attach itself to the cluster that minimises the free energy at that 
temperature and density. There may be more than one unique choice. The resulting change 
in energy and entropy of the system may be decomposed as the sum of two parts: a kinetic 
parti AE kin, A Skin), '^o with the cluster's motion in the environment; and an internal 
pa.Tt{AEin, ASin), to do with the internal motion of the particles constituting the cluster. 
If the volume is large, one may assume that the clusters form an almost ideal gas. In this 
case the average kinetic energy of a cluster of size k is (3/2)T. It does not depend on k and 
thus AEkin = 0. The change in internal energy AEin may be estimated simply as, 

AEin = -W + — + a{T) [{k + 1)2/3 - 0^] . (7.2) 
Co ^ 

The change in internal entropy is given simply as ASin = 2T/eQ. The kinetic entropy of 
an ideal gas of Uk clusters is given as(see Ref. pT|]), 



Skin f^k 



3 3 27rm 3 3 

- logT + log V + - log + - log A; + - 



-lognfc!. (7.3) 
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In most cases in nuclear fragmentation, Uk lies between and 1 (see Fig. (|TT])). Thus we 
may ignore the n^! term. Thus we get the total change in entropy for the addition of one 
particle to a cluster of size k as, 

A^ = 2T/eo + ^log(l + lA). (7.4) 

As a result, the total change in free energy and hence fi, is given as 

fx = AF = -W - — + a{T) Uk + 1)2/3 - P/^] - |t log (l + l/k) . (7.5) 



eo 



2 



We note that fi becomes progressively more negative with rising k. Thus the added 
particle prefers to attach to large clusters. In a fragmenting system under the boiling tem- 



perature, such a large fragment exists, of about half the size of the system [T^. The new 
particle thus preferentially attaches to this cluster. To illustrate this point more quantita- 
tively, we calculate this fi for a system with A = 240 and a.V = 4Vo. We assume the largest 



cluster is of size A/ 2. No doubt, this size falls gradually with rising temperature |10|, with 



the fall becoming rapid near the boiling temperature. In Fig. (12) we plot the value of 
this fi (solid line) assuming that the largest cluster remains of the same size throughout the 
temperature range. Above the boiling point the system exists mostly as small clusters, here 
we assume that the added particle attaches itself to a cluster of size k = 2. This fi is plotted 
as the dot-dashed line in Fig. (|T^. 



We are interested in obtaining an approximate expression for (n^) underneath the boiling 
temperature. Using the expression for fi as derived in Eq. ( [7.5[) , we obtain the expression 
for (rifc) as 

X exp [{(T(r)((A;^„. + if' - ^J^l) - ^Tlog(l + l/kma.)}k/T - a{T)e/' /t] 

For most systems, in general, the behaviour of kmax with T and V is difficult to estimate. 
However, we note that the above equation is precisely of the form of Eq. ( |7.1| ). On fitting the 
data points obtained from Eq. ( p.5|) we obtain the fit parameters as C* = 2.73, a = 0.36, 7 = 
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—3.11 (note that, as in the droplet model, only range of k between 10 to 40, is fitted). A 
plot of the fit to the values of ri^ obtained from a system with A = 240, T = 5MeV and 
V = 4Vo is shown in Fig. (0). Here the entire region from k = 2 to 240 is plotted. Note 
that both fits coincide extremely well in the region of /c = 10 to 40. 



VIII. SUMMARY AND DISCUSSION 

This investigation started out with the following question: suppose "experimental data" 
are given by the predictions of a theoretical model which, we know, does not conform to 
Eq.(l.l) exactly. Could we still describe the "data" approximately with the formula? And, 
if so, what are the significances of the parameters r, a and Tc? Subsequently we tried a fit 
with the droplet model and again found that a very adequate fit can be found. 

We find that the scaling law is still approximately obeyed. While we can not attach 
much significance to the extracted values of r or a (they are different from those given by 
the percolation or Ising model, but not terribly so), Tc seems to be a genuinely physical 
parameter, namely, it reflects the first order phase transition temperature. Coulomb effects 
tend to somewhat spoil even this correspondence. 

Many different fits can be obtained because by necessity the mass number of the com- 
posites is limited on the lower as well as on the higher side. We have no control over that 
since the dissociating systems are extremely finite. This apparently makes even deciding on 
the order of phase transition very difficult. Fits to data from calculations on percolation 
model have been made in the past . Since this is a model of continuous phase transition 
one necessarily concludes that the phase transition is continuous. It is more appropriate to 
consider the Lattice Gas Model [Q] instead which does encompass the Percolation model as a 



subset [^]. In the Lattice Gas model, if it is assumed that that the freeze-out density is less 
than half the normal density, then the fit of Eq. (1.1) would imply that Tc is indeed a first 
order phase transition temperature. If, on the other hand, one assumed that the freeze-out 
density is higher than half the normal density (we consider this highly unlikely) , then a fit to 
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Eq. (1.1) would not imply any usual thermodynamic phase transition pO[. If one depends 
on theories to decide on what order of phase transition to expect, one is driven towards 
expecting a first order phase transition. To date all the models which used a Hamiltonian 
0-0,|TO[ suggest a first order phase transition. 
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FIGURES 

FIG. 1. r, C^jA (left-hand scale) and (right-hand scale) plotted against temperature in the 
model of one kind of particles. The different panels are for different choices of A and V . 

FIG. 2. In < Ha > vs.lna. The solid line is the best fit to In < Ua > at each Tmax{a) presented 
by diamonds. The dotted line joining stars represents the distribution at Tc. The diferent panels 
are for different choices of A and V. 

FIG. 3. The scaling behavior in the mass range (10 < a < 40) 

FIG. 4. Same as Fig. 1, but in a model of two kinds of particles. Coulomb interaction has been 
switched off. The upper and lower panels are with respect to charge and mass distributions. 

FIG. 5. Upper panels: In < Uz > vs.lnz, and lower panels: In < na > vs.lna; in the model of 
two kinds of particles with Coulomb interaction switched off. 

FIG. 6. Scaling behavior in the charge range {7 < z < 17) (upper panels) and mass range 
(10 < a < 40) (lower panels) in two kinds of particles model without Coulomb interaction. 

FIG. 7. Tz and Cy/A as a function of temperature in a model of two kinds of particles, with 
Coulomb interaction included. 

FIG. 8. Same as upper panels of Fig. 5, but with Coulomb interaction. 

FIG. 9. Same as upper panels of Fig. 6, but with Coulomb interaction. 

FIG. 10. The parameters of droplet model a and 7 as a function of temperature for a system 
with A = 240 and V = 4Vo. The right pannels show the fit of the droplet model to the yields 
obtained in the model of one kind of particles described in section II. On the graph one can not 
distinguish between fitted points and the actual points from canonical calculations. 
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FIG. 11. < na> vs. a on a logarithmic plot. The right pannel expands the region 10 < a < 40. 

FIG. 12. n vs. T for a system with A = 240 and V = 4Vb. Data points represent results from 
a canonical calculation (see section II). Solid line represents /x for addition to the largest cluster, 
dot-dashed line is jj, for addition to a small cluster (see section VII for details). 

FIG. 13. Fits to < n<i > vs. a from two different models: open circles are from an exact 
canonical calculation; the solid line represents the fit by a droplet model; the dotted line represents 
the fit from the grand canonical approximation. 
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